Surface elastic waves in granular media under gravity 
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Due to the non-linearity of Hertzian contacts, the speed of sound c in granular matter is expected 
to increase with pressure as P 1 ^ 6 . A static layer of grains under gravity is thus stratified so that bulk 
waves are refracted toward the surface. We investigate wave propagation in the framework of an 
elastic description taking into account the main features of granular matter: non-linearity between 
stress and strain and existence of a yielding transition. We show in this context the existence of 
waves localised at the free surface - like Rayleigh waves - but with a wave guide effect related to 
the non-linear Hertz contact. The dispersion relationship shows a discrete number of modes which 
correspond to modes localized in the sagittal plane but also to transverse modes. The propagation 
speed of these waves is finally compared to recent measurements performed in the field, at the 
surface of a sand dune. 

PACS numbers: 



INTRODUCTION 

Desert animals like scorpions, use sand-born surface 
sound waves up to a distance of half a meter in order to 
localize their preyQ]. On each of their legs they possess 
slit sensilla receptors and associated neural connections, 
such as to detect small phase lags and orient their killing 
jump very efficiently. Actually to our knowledge, the 
first experiments on wave propagations over sandy free 
surfaces were conducted in this context. Biologists like 
Ph.Brownell and his collaborators |l| have identified sur- 
face propagation of sound and measured a rather low 
speed (around 50m /s in the 100-500 Hz range) com- 
patible with the animal biological capacities for signal 
processing^. From the physics point of view, surface 
waves in an elastic material (called Raleigh waves) are 
traveling at speeds slightly smaller than the bulk shear 
waves |3|- Nevertheless, in the context of granular mat- 
ter, the very existence of such waves and their ability to 
propagate over large distances is quite problematic. 

Granular matter, as a collection of stiff grains under 
moderate confining pressure, bears inter-granular con- 
tact surfaces of a scale much smaller than a typical grain 
size. In fact, the contact surfaces and the contact ori- 
entations are likely to be modified either reversibly of 
irreversibly under the variation of inter-granular forces. 
This feature is at the origin of many inherent difficulties 
when one wishes to determine the macroscopic constitu- 
tive properties. Moreover, the packing topology which 
determines the contact geometry is usually of a strongly 
disorderered nature |4J and average quantitities like de- 
formation fields are quite subtle to define 'Sj. Thus, al- 
ready at the most simple level of description, involving 
local non-linear elastic relations (the classical Hertz force 
problem 6])), a mean-field approach which identifies the 
local granular displacements with macroscopic deforma- 
tions, is failing quantitativelyjl H E3, O El E3 • 
In most cases, tests were made using sound wave bulk 
propagation^ OH, III 03 

under rather large confin- 




FIG. 1: Theoretical set-up. We consider elastic waves prop- 
agating along the x direction in a cell full-filled of granular 
matter. Gravity is along the z axis. The y axis is the direction 
transverse to the propagation. 

ing pressures P and results show a propagation veloc- 
ity c oc P 1 / 4 instead of the standard mean-field predic- 
tion c oc P 1 / 6 . Actually Makse et al.[12| have shown 
clearly using numerical simulations that if one takes into 
account the effective increase of the number of contacts 
with pressure, the agreement is bettered (especially for 
compression waves). Nevertheless, an essential discrep- 
ancy still lies in the assessment of the shearing stiffness 
merging from local tangential contact forces |l2|, LL3J ■ Note 
also, that from an experimental point of view, the ex- 
act origin of the discrepancy is not totally clair dial ). 
Other features such as a angular shape contacts jlaor 
the existence of a soft layer surrounding the grains |l 71] . 
can modify the propagation velocities in directions ob- 
served experimentally. 

In the limit of vanishing confining pressure very 
few results on sound wave propagation exist either 
experimentally Q, or theoretically |2jj. There are a 
priori several difficulties due to the unilaterality of gran- 
ular contacts which implies that the elastic constants are 
bound to vanish at zero pressure. Therefore, the elastic 
problem will go to a singular limit. A second order- like 
phase transition was recently identified in this limit, for 
frictionless grains around the J-point of jamming |2lj also 
displaying an anomalous density of low frequency vibra- 
tional modes j^. The study was extended to packing 
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with frictional contacts 23]. In practice, if one excludes 
microgravity experiments under controlled confinement 
(yet to be performed!), this peculiar situation would nat- 
urally exists close to a free surface. A simple (mean- 
field) arguments, shows that under these conditions, bulk 
waves cannot propagate parallel to the surface. This is 
due to the increase of material stiffness with depth which 
should redirect the wave front toward the free surface (the 
mirage effect). Experiments conducted by Liu et al.tlgj 
have shown an extreme sensitivity of the acoustic signal 
to minute local reorganisations, possibly due the speckle 
effect [1 4] , a dominant feature for probe sizes on the order 
of one grain. 

In this paper, we address the issue of surface wave 
propagation in a theoretical framework of continuous 
non-linear elasticity. The model we use was recently in- 
troduced by Jiang and Liu [24| to describe granular con- 
stitutive properties. In the first part of this paper, we 
give a short state of the art concerning the elastic descrip- 
tion of granular media in the static phase and present 
the phenomenology introduced by Jiang and Liu. In a 
second part, we derive the equations of propagation for 
an isotropic compression and and uni-axial load. In the 
third part, we analyse the propagation of Ray leigh- like 
surface waves under gravity and compare the results to 
experimental measurements made in the field, on a dune 
slope. 

THE JIANG-LIU MODEL 



many qualitative features observed experimentally, such 
as the existence of a Coulomb-like failure or stress in- 
duced anisotropy. Efforts have been made by the authors 
to compare quantitatively their model to the output of 
several experimental measurements, such as systematic 
triaxial tests, the response to a local load and static equi- 
librium in a column. The agreement was noticeable, es- 
pecially in view of the minimal amount of free parameters 
in the model. The possibility to follow the elastic behav- 
ior up to the limit of failure is also a promising feature of 
the model in the context of slope stability monitored by 
sound waves. Note that this feature is quite original if 
one compares to the standard Boussinesq non-linear elas- 
ticity framework [ij|(see ref.j^). Of course, issues like 
irreversible deformation fields (plasticity) are still ques- 
tionable in the framework of this model; neverthelessc it 
provides a well defined starting point for a complete anal- 
yses and modelling of elastic vibrations. Moreover, this 
approach can be generalized easily with a power law not 
representing necessarily the Hertzian interaction. Thus, 
one could in principle, take into account the existence of 
other types of non-linear contact force-laws E1II3. In 
the present paper, we will limit ourselves to Hertzian in- 
teractions and we will derive sound propagation around 
the simplest possible reference states. 

The derivation of the stress tensor yields : 

<7« = EVS (BSS^ - 2Au% + ^p^j (4) 



We introduce the displacement field pU and the stress 
tensor oij. With our conventions, the dynamical equa- 
tion reads: 



pU = -Vcr + pg 



(1) 



The stress tensor cry is related to the strain tensor ity — 

1 (dUt | dUj N 

2 \dxj ' dxi / 

Recently, Jiang et Liu have proposed a phenomenological 
formulation for the elasticity of granular states based on 
a simple energetic formulation : 



F eI = ^ES 1 / 2 (Bu u S 2 + 2i 



(2) 



where E is the material young's modulus, A and B are 
two dimensionless numbers. The volumic compression is 



-Tr{uij) 



(3) 



U Q tj = Ui 



^Sij is the traceless strain tensor and 



its modulus squared. Interestingly, this elegant and com- 
pact formulation of elastic energy is able to reproduce 



Mean field determination of the compression 
modulus 



Since Mindlin in the 50's, several works, focusing 
on sound wave bulk propagation in granular materials, 
have derived from local Herzian interactions, the effec- 
tive elastic moduli 0, IE El Hof ; for a recent and clear 
review, see Makse et al..l3j and references inside. We 
just recall here the mean-field results for a packing of 
spheres of compacity $ and an average coordination num- 
ber Z. First let us write the interparticular forces be- 
tween two spheres of radii R\ and R2 and half overlap 
£ = 1/2(R 1 +R 2 - I z\ - ^ 2 |), R = 2R 1 R 2 /(R 1 + R 2 ) 
(if R\ = R 2 , R is the sphere radius). The normal force 
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where p 9 is the shear modulus and v g the Poisson ratio. If 
half the tangential separation between the sphere center 
is As, the tangential force is then : 



8p s 



-R 



1/2 



As 
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Finally, one may obtain from mean-field granular 
displacement, the elastic constant values, i.e. the bulk 
modulus: 



1/3 



37r(l-2/g) V 2 Ms 

and the shear modulus: 



1 — Vg 2 — Vg 5lT \ 2(lg 
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with two limits (i) no sliding friction between the grains 
(take e = 0) or (ii) infinite friction between the grains 
(take s = 1). For the macroscopic Jiang-Liu model in 
the isotropic case, the strain and stress tensors read : 

mj = a ij = Phi P = EB5^ /2 (5) 

3/2 

For an isotropic compression, we have: P = EB5 Q — 
Kmf^o, therefore: 



K MF = EE 



1/3 



Consequently 



B 



2 3 / 2 3tt(1 - i/3) 



(6) 



Determination of the ratio B/A 

Mean field 

Now we consider the Jiang-Liu model with x-z shear 
under isotropic compression P : 

1/2 



<?xz = -2AE5 ' = -2(j, MF u x 



and 
Thus, 

Therefore : 
A 



P = EES, 



3/2 



[IMF = EA 
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2 3 / 2 5tt(1 + ^) \l-v g 2- 



1 3e 



(7) 



If one takes $ ~ 0.6 and Z ~ 6. The bulk shear modu- 
lus for silica oxide is /i g ~ 30GPa and the Poisson ratio is 
v g ~ 0.2. This gives : £5 ~ lOGPa and EA ~ 6.5GPa 



(for e = and £M ~ 9GPa for e = 1. These values will 
be of importance when the surface sound wave velocities, 
found experimentally for sand, will be discussed at the 
end of the paper. 

From the expressions above, we get: 

(8) 



A 



(1 + 3^) 
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For e = 0, we get B/A = 5/3. For e = 1 and v g 
we get B/A = 5/7. 



0.2, 



Energetic argument 

According to Jiang and Liu, we have a material insta- 
bility corresponding to a Coulomb yield, for a criterion 
based on the free energy landscape convexity (a thermo- 
dynamic stability criterion) |24j . 



tan# = 



(9) 



Thus B/A =2/ tan 2 9 C ~ 6 We then notice a large de- 
screpancy for the B/A ratio between the mean-Field so- 
lution and the empirical result obtained from direct mea- 
surement of the sand-pile slope. 



SOUND PROPAGATION UNDER ISOTROPIC 
COMPRESSION 

Longitudinal waves 

Now, let us consider the propagation of longitudi- 
nal waves along the x axis. As the system is homoge- 
neous, the modes are simply Fourier modes of the form: 
Tj e t(kx-ut) _ denote by the disturbance to the 

strain field and ul 



-o • 



ikU 




its traceless counterpart 




ikU 
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The disturbance of the volumic compression is then: 6 = 
—ikU and that of the modulus u 2 s is null. The stress 
associated to the sound wave is: 



= Ey 



\b55h 



2Au l 



so that the equation of motion finally reads: 
- pu 2 U = -k 2 Eyf^ 



u 



The speed of longitudinal acoustic waves is finally: 

1/6 



T|| 



(10) 



(11) 
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with 7|| = ±A) 1/2 B- 1/6 . We thus recover the 

scaling of the speed of sound c as P 1 ' 6 , but with an extra 
dependence on the coefficients A and B which themselves 
depends on the mean number of contacts and thus on the 
pressure. 

Transverse waves 

Similarly, we consider the propagation of transverse 
waves along the x axis. The modes are still Fourier modes 
of the form: Ve i( - kx - ut h 

/ 1/2 0\ 
ufj = ikV 1/2 
\ / 

The disturbance of the volumic compression 5 and of the 
modulus u 2 are both null. The disturbance to the stress 
reduces to: 

o7j = -2Ey/6oAuij° 
and the equation of motion finally reads: 

- puj 2 V = -k 2 E^j%AV (12) 
The speed of longitudinal acoustic waves is finally: 




with 7_l = A 1/2 B- 1 / 6 . 

We thus get a prediction for the relation between the 
propagation speeds of longitudinal and transverse waves: 

Jl = (&L + lV /2 = ( J_ _ 1 V /2 

7± \2A 3J \2p 2 6J 



TRANSVERSE SURFACE WAVES UNDER 
GRAVITY 

Equilibrium of the system 

The question is now to examine the case of a large 
rectangular box filled with granular matter and submit- 
ted to a vertical gravity field (fig. Q). As mentioned in 
the introduction, due to the inhomogeneous pressure field 
P oc pgz, the system presents a stratification in the prop- 
agation speed c oc z 1 / 6 which should also lead to a refrac- 
tion toward the free surface. At first, it is clear that this 
mirage effect does not affect the propagation along the 
z direction. Consequently, due to reflections also occur- 
ring at the free surface, we can expect that the box will 
present a series of vibration resonances. 




FIG. 2: Shape of the transverse modes for the same frequency. 
The modes n > are cut 

The first step is to compute the elastic equilibrium 
state of the system under gravity. We first solve the 
equilibrium problem, starting from a strain tensor of the 
form: 

/ \ / 1 \ 

Uii = u% = 5 | 

V -So J \ -f / 

The modulus follows as u 2 = |(5g. Using the constitutive 
relation, we obtain the stress tensor: 

<7 t ,=E5 3 /2 [ (B-i) 

V (B+¥)J 

The pressure is defined from the trace of this stress ten- 
sor: P = ES 3 /2 (B+£) The expression for the vertical 
stress is replaced by: 

^=(B+^E5r = P9z 
so that the volumic compression reads: 

dW f pgz \ 1 

where Wq is the vertical displacement field. At the lower 
edge of the box z = H (fig.[TJ, we impose that the grains 
remain fixed: the boundary condition is Wo(H) = 0. 

Linearized equations 

Now, we investigate the existence of - propaga- 
tive - transverse modes localized close the surface(see 
ref.[2(|). We consider a transverse displacement <^(a;, £) = 
^ e i(kx~u)t) a |. ^ e sur f ace f granular bed. We define 
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n=4 n=5 i oo-i 




FIG. 3: Shape of resonant - not propagative - modes. 

the dimensionless function V of the dimensionless vari- 
able rj = kz: 



U = | V(kz) | c 
Again, we linearize the strain tensor: 



(15) 



,.v 



u$ = K\ f o £ 



£ 



and find no variation of the volumic compression (5 = 0. 
The trace of Uij is null so that u\j = Uij, and then 
u s 2 — 2uijUij vanishes. Finally we obtain : 

iV 

-2AE^5 ufj = -AEy/dokC | iV V 

V 

Again, from the equation of motion along y, we get a 
dispersion relation within a constant a: 



E 



1/3 



uj = aA 1/2 (B + 5A/3)- 1/6 i^-j g x ' 6 k 5 ' 6 (16) 
The shape of the mode is given by the equation: 

a 2 V = r)V 3 V - (n^V')' (17) 
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FIG. 4: Rescaled pulsation a as a function of kH for Rayleigh- 
Hertz waves (a) and transverse surface waves (b). c) Ratio 
as a function of, for transverse surface modes. The different 
branches collapse on a single curve. 



with the boundary conditions coming from the zero stress 
condition at the free surface and the definition of £: 
V(0) = 1, V"(0) = 0, V(kH) = 0. 

For infinite depth, the system presents a discrete num- 
ber of modes which may be associated to the stratifi- 
cation. The figure shows the variation of otoo with n 
together with a fit by (oto + San) 1 / 6 . We have computed 
systematically the ratio a/aoo as a function of kH. For 
small kH, we find that a decreases like (kH)~ 5 ' 6 . For kH 
much larger than 1, a tends to Qk,. The curves for the 
different modes collapse when kH is properly rescaled. 



Resolution 



Wave-guide cut-off: horizontal resonant modes 



We decompose the problem into two equations: 

v = v - 1/3 s 

S> = _ a 2 )V 



The height of the box is the only relevant length scale 
of the problem. We introduce accordingly the coordi- 
nate r\ = z/H. We consider modes of vibration whose 
disturbed displacement field is horizontal, of the form 
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U(z/H)e 



The disturbed strain field reads: 



and is traceless counterpart is: 



U'(r])/2H 

U'(rj)/2H 



and find no variation of the volumic compression S = 0. 

o — , ~ 2 
= Uij, and u s — 

-2AEy/tiEuTj- The 
disturbed vertical stress deduces as: 



The trace of u^j is null so that u ZJ 
vanishes. We end with a i3 = 



o7 x = -AE^/SoU'/H 
and leads to the equation of motion: 

(r, 1/3 U'(v))' + P 2 U( V )=0 
where (3 is the rescaled pulsation defined by: 



(18) 



PA 1 ' 2 ( B 



5A 
~3~ 



-1/6 



1/3 



g l/6 H -5/6 ( 1Q ) 



Remarkably, the resonant frequencies are found to scale 
as _ff -5 / 6 and not H~ x as for a non dispersive medium. 

The solution of the ordinary equation (|29fl involves the 
Bessel function of first kind J 2 /5 and is of the form: 



U(rj) = rj 2/3 J„ 



2/5 



5 



(20) 



The value of (3 is selected by the zero displacement con- 
dition at the bottom edge of the box, U(l) = 0, which 
simplifies into J2/s(6/3/5) = 0. There is thus a discrete 
number of modes, labelled by n, whose rescaled frequency 
is: 



6/5 



5nn 
~6~ 



(21) 



RAYLEIGH-HERTZ WAVES UNDER GRAVITY 
Linearized equations 

Now we investigate the existence of propagative modes 
analogous to the Rayleigh modes i.e. whose deformations 
are in the sagital plane (x, z). We consider the defor- 
mation of the free surface to be of the form ((x, t) = 
( e ! ' tl_1Jt '. The displacement field is defined trough two 
dimensionless functions U and W of the dimensionless 
variable r/ = kz: 



U = 



The disturbance to the strain tensor reads: 






f -u 
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1 i(U' + W) 
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/ _2U _ n i(U' + W) 

3 3 2 
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. i(U' + W) n 2W^_ , U 

\ 2 u 3^3 



The disturbance to the volumic compression is 6 — 
k((U — W). That of the modulus u 2 s reads: 



u s 2 = 2u°,u°j = 



2<5 



(U + 2W')k£ 



The disturbance to the stress tensor may be formally 
expressed as: 



'SaE 
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which gives, after simplification, the following expres- 
sions for its components: 

S£ = EVK [(i + f ) u + (# - f ) W] k( 
^ z = - Ey fc [(4 - f ) u + + f ) W>] k( 

^r z = ^7 x = -iAEVh* {u 1 + w) k 

For the sake of simplicity, we rescale the stress tensor, 
introducing the functions S xz (kz), S zz (kz), S xx (kz): 

^ = liEkf^ipgf^AiB + 5A/3)- 1/3 (S xz {kz) 
<T ZZ = {Ekf^ipgf^AiB + bA/Zy^tSUkz) 
^ = {Ekf^ipgf'^AiB + 5A/3)- 1/3 <:S xx (kz)(22) 

The previous expressions then simplify into: 



,1/3 



S Z M =-r? 1/3 
S X M =-r? 1/3 [U'{v) + W(r})] 



1 3B\ TT/ . (I 3B , 

1 3B\ TT , . /5 3B , n , 
2-2a) U ^ + \-2 + 2A^ W ^ 



(23) 



From the equations of motion, we get the dispersion 
relationship within a constant a: 



aA 1 ' 2 {B + hA/Z)- 1 ' 6 (jj <7 1/6 fc 5/6 (24) 



as well as a set of equations governing the shape of the 
modes: 



a 2 U — S xx + S xz 
a 2 W = -S xz + S' zz 



(25) 
(26) 
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FIG. 5: Dispersion relation of longitudinal and transverse 
waves. 



Resolution 

In order to solve the boundary condition problem, we 
rewrite the equations as a set of 4 - linear - ordinary 
differential equations: 

U' = -W-rj-^S.^ 
1 



w 

^xz 
$zz 



hA + ZB 

a 2 U — S xx 
o?W + S xz 
1 

5A + 3B 



(A -3B)U + 2Arj- 1/3 S z 



{3B - A) S zz + 2{A + 6B) rj^U 



with the boundary conditions: U(kH) = 0, W(0) = 1, 
W(kH) = 0, S xz (0) = and ,5* 2Z (0) = 0. The solution 
[U, W, S xz , S zz ](ri) is obtained by superposition of the 
solutions obtained, starting from the initial conditions 
[0, 1, 0, 0] and [1, 0, 0, 0]. The value of a is then tuned to 
get the cancellation of both U(kH) and W{kH) = 0. 
The central result of this part is that contrarily to 
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FIG. 6: a Asymptotic value of the rescaled frequency a in 
the limit of large kH , as a function of the mode number n. b 
Asymptotic value of the rescaled frequency /3 in the limit of 
small kH. 



the classical Rayleigh waves, the propagation of these 
" Rayleigh Hertz" surface waves presents a discrete num- 
ber of modes that may be associated to the stratification 
effect. The shape of the modes excited at a given pul- 
sation lo is displayed on fig. [7] For infinite depth i.e. in 
the limit of infinite kH, the rescaled pulsation a tends 
toward a constant aoo whose dependence on the mode 
number n is shown on fig. There is a general trend of 
ctoo to increase as n 1 / 6 but the values are not regularly 
spaced. 

For small kH, we find the usual wave-guide cut-off as 
lo tends to a constant. In terms of a, this corresponds to 
an asymptotic behavior of the form: 

a = (3(kH)- 5 / e (27) 

So, in the limit of wavelength large compared to the cell 
depth, the pulsation tends to the value 



u> 



PA 1/2 (B + 5A/3Y 



1/3 



V6ff-B/« (28) 



The value of j3 as a function of n is displayed on fig. [S] 
Although it is not a regularly increasing function of n, (3 
follows a general trend of being linearly dependent on n. 



Wave-guide cut-off: vertical resonant modes 

We now consider modes of vibration whose dis- 
turbed displacement field is vertical and of the form 
W(z/H)e lwt . Then the disturbed strain field reads: 
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FIG. 7: Shape of the Rayleigh-Hertz modes for a same fre- 
quency. The modes n > are cut 



from which we deduce the disturbed volumic compression 
S = —W'/H and the disturbed modulus 



u, 2 = 2v2*u? - 



4S W' 



Then the disturbed vertical stress is deduced as: 



2H 



(5A + 3B) W 



and leads to the same equation of motion as previously 
(V /3 W'(Tijj + f3 2 W(r)) = (29) 
but where (3 is now defined by: 

1/3 



uj = (3 



(30) 



EXPERIMENTAL MEASUREMENTS AND 
DISCUSSION 

In conclusion, we performed a derivation of sound-wave 
propagation modes in the context of the Jiang-Liu model 
of granular non-linear elasticity. First, the longitudinal 
and transverse sound velocities for and isotropic compres- 
sion were calculated. Morerover, we derived the propa- 
gating and non propagating modes for a granular packing 
with a free surface under a vertical gravity field. We show 
that the gravity field producing a stratification of the 
material stiffnes is responsible for a channeling effect of 
the acoustic waves. Hence, we find two surface propaga- 
tion modes in the transverse and saggital directions. The 
last one, we call Rayleigh-Hertz waves. Now we present 
the results of experimental measurements performed on a 
sand dune using localized surface excitation and already 
reported in reference |27|. Sinusoidal signals, recorded 
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FIG. 8: Dispersion relation of surface elastic waves on a dune 
slip face. Transducers distant by 5 cm (down triangle), 15 cm 
(up triangle), 25 cm (square) and 42 cm (circle). The solid 
lines corresponds to the saggittal Rayleigh-Herz modes (see 
equ. ) with a material stiffness adjusted to match the lowest 
n = propagating mode. 



on a tape, were played through an amplified loud-speaker 
that was kept at the surface of the avalanche slip face of 
a singing dune (top schematic of fig. |SJ. 

As the mechanical transmission of the sound gener- 
ator with the sand bed is non-linear, the signal is dis- 
torted and presents harmonics. This imperfection was 
used to derive several points of the dispersion relation 
from each recording. From the phase between the Fourier 
components of the signals of two accelerometers aligned 
with the loud-speaker, the wavelengths A associated to 
each frequency / were obtained. In order to unwrap the 
phase, the amplitude of the signals was slightly modu- 
lated at low frequency. On fig. we present the col- 
lection of all measurements that establish the disper- 
sion relation. Actually, when replaced in the context 
of Rayleigh-Herz surface waves, the present framework 
sheds a new light on the apparent large scatter of the 
data. Indeed, we can superpose on the experimental 
results, the dispersion relations found for the sagittal 
vibrations. On the figure, the solid lines correspond 
to the dispersion relations : / = c„fc 5 / 6 , with c„ = 
a aD {n)A 1 / 2 (B + bA/'i)~ 1 / Q (E/ P ) 1 / :i g 1 / & corresponding to 
the n th mode prefactor derived for an infinite depth (see 
fig. EH for aoo{n) values). By best fit, the Co value was 
adjusted to match the graviest mode and we obtained 
Co ~ 50Z7 SI. Note that the agreement for n — 0, 1 and 
2 is quite good. If we compare this result to the values 
derived from the mean-field estimation EB = lOGPa 
and the B/A = 6 derived from the sand-pile slope, we 
get c ~ 105USI. Which means that the actual surface 
waves speeds measured in these conditions are still slow 
by a factor 2 when evaluated from the simple arguments 
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presented above (or alternatively would correspond to an 
overestimation of A by a factor 4) . Thus, in spite of this 
semi-quantitative agreement, there are still many ques- 
tion concerning the direct evidence of the propagation 
modes derived in this context. This is an ongoing work 
in our laboratory and we leave these questions for a fu- 
ture report. 
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